Experimental investigation of pair dispersion with small initial separation in 

convective turbulent flows 
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We report an experimental investigation of pair dispersions in turbulent thermal convection with 
initial separation ro ranging from sub-Kolmogorov scale to scales in the inertial range. In the 
dissipative range of scales we observed for the first time in experiment the exponential growth of 
the separation between a pair of particles predicted by Batchelor and obtained a Batchelor constant 
0.23 ±0.07. For large ro, it is found that, for almost all time range, both the mean-square separation 
and distance neighbor function exhibit the forms predicted by Batchelor, whereas the two quantities 
agree with Richardson's predictions for small ro. Moreover, the measured value of the Richardson 
constant g = 0.10 ± 0.07, which is smaller than those found in other turbulence systems. We also 
demonstrate the crossover of the mean-square separation from the exponential to the Batchelor 
regimes in both temporal and spatial scales. 

PACS numbers: 47.27.-i, 44.25. +f, 47.55.pb, 47.27.tb 
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Turbulent relative dispersion of a pair of particles is of 
central importance to a wide range of natural processes 
such as pollutant spreading in the atmosphere [l[ and 
mixing in oceans. The concept was first introduced by 
Richardson, who attempted to explain the large observed 
value of turbulent diffusivity in the atmosphere. He in- 
troduced a quantity named distance neighbor function 
(DNF) and the diffusion equation to describe the evo- 
lution of DNF 0]. With Kolmogorov's scaling theory, 
Obukhov refined Richardson's prediction and found that 
(r 2 ) = get 3 , with r the pair separation, e the mean kinetic 
energy dissipation rate, and g is a dimensionless constant 
called the Richardson constant Q. Batchelor recog- 
nizing that over short time the initial separation ro be- 
tween the pair of particles would be important, obtained 
(|r(t)-r | 2 ) = f(r Q )t 2 for t v « t « t , where t v = (is/e) 1/2 
is the Kolmogorov time scale and to = (rg/e) 1 / 3 is a char- 
acteristic time below which the initial separation is im- 
portant. The Richardson-Obukhov scaling is now sup- 
posed to hold for to « t « tr,, where ti is the integral 
time scale. In the above f(r ) = [D LL {r ) + 2D NN (r )], 
with Dll and Dnn being the 2nd-order longitudinal 
and transverse Eularian structure functions (ESF) re- 
spectively. 

An important regime of pair dispersion is the very early 
stage of separation in which relevant spatial scales are 
within the dissipative subrange [1, [f| . Turbulent disper- 
sion and mixing in this regime is closely related to the 
reaction rate for fast reacting scalars, such as in combus- 
tions. Batchelor was the first to argue that the growth 
rate of pair separation should be proportional to the sep- 
aration distance r itself in this regime Q , which leads to 
an exponential growth (r 2 ) ~ r^exp^t) when both the 
initial and final particle separations are within the dissi- 
pative range, r 2 , « {r 2 (t)) « rj 2 . Here the growth rate 
£ = and B is called the Batchelor constant. How- 

ever, in most previous experimental studies, the initial 
separation is in the inertial subrange @, As a result, 
the exponential growth regime has never been observed in 



experiments and the Batchelor constant has never been 
measured. 

Recently, we have shown that the Lagrangian particle 
tracking vclocimetry (PTV) can be applied to thermally- 
driven turbulent flows and have obtained particle pairs 
with separations smaller than r\ (Io| . One advantage of 
our system is that the range of its parameters is such 
that both the dissipative and inertial subranges can be 
easily accessed in the experiment. In fact, we have accu- 
rately determined the energy dissipation rate e from the 
measured dissipative range ESFs using PTV [ll| . In this 
respect, turbulent thermal convection provides a good 
platform for studying properties of particle dispersions 
in both the dissipative and inertial subrange in a sin- 
gle experiment. From a practical point of view, studying 
two-particle dispersion in turbulent thermal convection is 
important in understanding the motion of passive scalars 
in the atmosphere and oceans in which buoyancy is a rele- 
vant driving force and is absent in most previous studies. 

In this Letter, we report measurements of the mean- 
square separation of a pair of particles (r 2 (t)) (and also 
(|r(i) -ro| 2 }) in buoyancy-driven turbulent thermal con- 
vection. Lagrangian particle tracking velocimetry (PTV) 
was used in the experiments, which were carried out in a 
cylindrical cell with water as working fluid. The measure- 
ments were made in the cell center with ro ranging from 
dissipative to inertial range of scales. The height and di- 
ameter of the cell both equal to 19.2 cm, so the aspect 
ratio is one. The experiments were conducted at fixed 
Prandtl number Pr = ^/k=4.4 with various Rayleigh 
number Ra = agATH 3 /vn (from 2.9 x 10 9 to 1.3 x 10 10 ), 
here g is the gravitational acceleration, AT the temper- 
ature difference across the fluid layer, and a, v and k, 
respectively, the thermal expansion coefficient, kinematic 
viscosity and the thermal diffusivity of the working fluid. 
To compare the results with other turbulence systems, 
the micro-scale Reynolds numbers R\ are determined by 
using R\ = y/l5u' 4 /ev with u' the root mean square ve- 
locity and e the energy dissipation rate in the cell center 
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fioj . The tracking volume [SV -(5 cm) 3 ] in the center of 
the cell was illuminated by a laser beam, and the scat- 
tered light from the seeding particles (diameter d p = 50 
/mi polyamid, density p = 1.03 g/cm 3 ) were captured 
by three cameras simultaneously. The Stokes number 
St = t p /t v ranges from 10~ 4 to 10~ 3 , with t p being the 
time scale of the Stokes viscous drag due to interaction 
between particle and fluid. The number is much less 
than one, indicating that the particles would be safely re- 
garded as tracers. As Kolmogorov time scale is 0.3~0.5 s 
in the parameter range of our experiment and the camera 
frame rate is 50 or 100 fps depending on the Ra, the tem- 
poral resolution is sufficient to resolve dissipative range 
properties. The error of the particle position after cali- 
bration is ~ 8 /im. However, due to the finite particle size 
and the diffraction effect, the minimum resolvable sepa- 
ration between a pair of particles is typically between 100 
/xm to 200 /xm, which is less than r\ ^ 0.5 mm in the ex- 
periment. Thus we are able to determine pair separations 
with initial separation ro smaller than the Kolmogorov 
length scale. In practice, we binned all pairs of particles 
with initial separation in the range [ro-<5r , tq+Stq] when 
computing statistics. For ro < 20?7, we take 5r « 77/5 and 
for r > 207/, Sr ra 77. Other details of the setup and 
calibration have been described elsewhere [l(| ■ 




FIG. 1: (Color online) The mean square separation (r 2 ) as 
a function of time for four different Rayleigh numbers with 
initial separation smaller or close to 77 in a semi- log plot. From 
top to bottom, cyan squares (Ra = 1.0 x 10 10 , R\ = 67, ro = 
77); red circles (Ra = 1.3 x 10 10 , -Ra = 84, ro = 0.9r/); green 
triangles (Ra = 6.1 x 10 9 , -Ra = 53, ro = 0.7r/); pink diamonds 
(Ra = 2.9 x 10 9 , R A = 35, r = 0.4r,). The dashed lines show 
the exponential fit to the respective data for t < t v . 



Figure [T] plots in semi-log scale four mean square sep- 
arations all with ro within the dissipative scale. It is 
seen that for t < t v and r < 1.5r/ all curves with differ- 
ent Ra (-Ra) do follow an exponential growth as shown 
by the respective dashed lines. This initial exponential 



growth is observed for all nine values of Ra measured. 
From exponential fits we find no systematic Rayleigh 
(Reynolds) number dependency for the Batchclor con- 
stant, which leads to an average Batchelor constant B = 
0.23±0.07. This number was first estimated by Batchelor 
and Townsend based on the assumption that the dissipa- 
tive separation is mainly due to stretching from the ve- 
locity gradient, which gives a range B = 0.35 ~ 0.41 [i"2| . 
However, it was later argued that the original estimate 
is too large due to the lack of persistence of the rate- 
of-strain tensor and the role of vorticity [l3| and it only 
serves as an upper limit. There are several simulation 
[l3| and model [l4[ studies that attempt to estimate this 
constant and they give B = 0.093 ~ 0.13. It is seen that 
our value of Batchclor constant is larger than previous 
findings but still smaller than the upper limit proposed 
by Batchelor & Townsend [l2|]. To our knowledge, the 
present result is the first experimental confirmation of 
the exponential regime and the direct measurement of 
the Batchelor constant. 
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FIG. 2: (Color online) The normalized mean square separa- 
tion ((r - ro) 2 )/r/ 2 as a function of time for different initial 
separations. From bottom to top, ro=0.9r/, 1.3r/, 1.7?7, 2.2r/, 
4.377, 6.57/, 8.777, IO.877, I3.O77, 15.1t/. The dotted line and 
dashed line indicate the Batchelor and Richardson regimes 
respectively. 



Figure [2] shows the temporal evolution of (\r(t) -ro| 2 } 
for different ro, which increases, from bottom to top, 
from dissipative range of scales to inertial range of scales. 
The dotted and dashed lines in the figure show Batchc- 
lor and Richardson scalings respectively In PTV, the 
number of velocity pairs varies for different spatial sep- 
arations, and pairs with the small separations generally 
have lower probability of being measured than those in 
some intermediate range of scales. Therefore, the one 
with smallest ro has lowest number of particle pairs for 
statistics. In our experiment, there are 10 pairs of par- 
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tides for ro = 60. 677, but only 10 3 for ro = 0.9?/. In 
Fig. [5] we show only the statistical errors for r = O.917, 
as the uncertainty for this data is the biggest. One may 
note that even the largest error bar is within the symbol. 
For large ro, there is a power law regime in a range of 
time scales extending from 0.1 to 6 t v , whose exponent 
is very close to the one predicted by Batchelor [J], i.e. 
(|r(i) -ro| 2 ) ~ t 2 (shown as the dotted line in the figure). 
For ?'o = 0.9?7, there is no single power law for the entire 
range. For t > 2r„ the behavior could be well described by 
the Richardson-Obukhov's law for two-particle diffusion, 
i.e. (|r(£) -i"o| 2 ) = get 3 . This can be seen more clearly in 
the inset of Fig. [3] as the plateau for the square separa- 
tion compensated by Richardson-Obukhov scaling. Here 
it is seen that two data sets with the smallest initial sep- 
arations (ro = 0.9 and 1 .377) reach the plateau, whereas 
the others show similar trend towards the Richardson- 
Obukhov scaling but lack sufficient time to develop. Note 
that the time scale in the inset of Fig. [3] is normalized by 
a characteristic time to = (ro/e) 1 ^ 3 , below which the ini- 
tial separation is important. From the red solid line we 
obtain the Richardson constant g = 0.10 ± 0.007. Pre- 
vious studies show that there is a large uncertainty on 
the value of g. For non-buoyancy driven turbulent flows, 
more recent experimental and numerical studies 0, f]~5l — 
HH suggest g « 0.5. For thermal convection, a numerical 
study found that the value of g equals to 0.16 and this 
smaller value was attributed to the correlated pair mo- 
tion in thermal plumes [l|l]. As the plumes' motions are 
predominantly in the vertical direction, this would imply 
that pair dispersions should behave differently in differ- 
ent direction. However, by studying pair dispersion in 
the vertical and lateral directions separately we find the 
dispersion properties to be isotropic with the Richardson 
constant nearly the same in the vertical and the two lat- 
eral directions, i.e. each being 0.03. This suggests that 
pair dispersion in all directions are affected by some cor- 
related motions. It is noted that flow in the cell center 
is affected by the large-scale circulation. This coherent 
motion is also azimuthally rotating, which may induce 
correlated motions in different directions. 

The DNF represents the spherically averaged PDF 
for pairs of particles with separation r at time t, 
i.e. p(r,t). Richardson first suggested that rela- 
tive dispersion can be modeled by a diffusion equa- 
tion for the DNF. For the isotropic flow, the dif- 
fusion equation can be expressed as dp(r,t)/dt = 
(l/r 2 )d[r 2 K(r,t)dp(r 7 t)/r]/dr with K(r,t) being the 
diffusion constant. Richardson proposed that K(r, t) = 
koe 1 / 3 ^/ 3 based on the experimental measurements 
in the atmosphere, and found that pji([r,t) 
v / 143/2x429/70(7r(r 2 ))- 3 / 2 exp[-(1287r 2 /8(r 2 )) 1 / 3 ]. As- 
suming K(r,t) ~ t 2 , Batchelor found another so- 
lution to the diffusion equation, i.e. PB(f, t) 
(2n(r 2 )/3)- 3,2 exp[-3r 2 /2(r 2 )]. The two solutions are 
shown in Fig. [3] as black solid line (Richardson's pre- 
diction) and red dashed line (Batchclor's prediction). In 
both solutions, the separation between two particles were 




FIG. 3: (Color online) The distance neighbour functions for 
different initial separations at Ra = 1.3 x 10 10 (R\ = 84). The 
red curved dashed line is Batchelors predicted PDF, while the 
black straight line is Richardsons. The open symbols show the 
experimental results for initial separation ro = 0.9r/ with time 
t ranging from 3. It,, to 5.5r, ? . The closed symbols show the 
results for initial separation ro = 52-q with time t ranging from 
5.5t v to 8.1T,,. Inset: The mean-square separation compen- 
sated by t 3 with time normalized by to- The red solid line 
gives the Richardson constant g = 0.01. From bottom to top, 
the curves represent different initial separations same as those 
in Fig. [2] The thickness of each line at different times show 
the uncertainties for the data. 



assumed to be zero at the very beginning. Experimen- 
tally, however, even one could resolve sub-Kolmogorov 
scale, the initial separation would be much larger than 0. 
One way to solve this problem is subtracting all particle 
separations with their initial value Ar = r-ro, and replac- 
ing p(r,t) with p(Ar,t) Q. In Fig. [31 the open symbols 
all have ro = 0.9rj and closed ones ro = 52r/. There are 
five data sets at different times for each initial separa- 
tion. The time are chosen to fall into the time range 
where the particle separations increase as t 2 (ro = 52r/) 
and t 3 (ro = 0.9r/) scalings respectively. It is clear that 
the DNF results agree with Richardson's prediction for 
small ro and agree with Batchclor's prediction for large 
ro- 

To take a closer look at the Batchelor regime, we 
show in Fig. 3] the mean square pair separation com- 
pensated by f(ro)t 2 . Note that because of initial sep- 
aration ?*o in our experiments varies continuously from 
the dissipative range to inertial range, we use the full 
function /(r ) = Dll^o) + 2£>ATAr(ro) for the coefficient 
of the Batchelor scaling, instead of its dissipative range 
(|ro/r/ 2 ) or inertial range (jCr^ 3 /j) 2 ' 3 ) scalings as in 
some previous studies. It is seen from the figure that 
curves for all initial separations and for t from r T; to 3t v 
collapse onto one horizontal line with the height very 
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close to unity. In the above the values of the ESFs /(ro) 
were independently obtained from the measured parti- 
cle trajectories which are also shown as the solid 
blue line in the inset of the figure. The values of /(ro) 
can also be obtained as the plateau heights of the com- 
pensated plots (|r(t) - ro\ 2 )/t 2 for various values of ro 
(not shown here), which are shown as the circles in the 
inset. It is seen that there is an excellent agreement be- 
tween the values of the /(r ) obtained from the measured 
mean square pair dispersion and the Bachelor relation 
and those obtained directly from ESFs. Also shown in 
the inset are the K41 predictions for the dissipative (red 
line), inertial range (green line), and the large ro limit 
6Rx/VTE (red dashed line) of f(r ) 0. The excellent 
collapse between circles and three solid lines indicates 
that the dispersion in the intermediate time domain is 
mainly controlled by the initial velocity difference be- 
tween two particles with their separation extending from 
dissipative to inertial ranges. In Fig. 01a) it is seen that 
the compensated (\r(t) - tq\ 2 ) for ro ~ f] and t < T n in- 
creases systematically as ro decreases. The reason is as 
follows. For very small t and ro < 77 (so the velocity dif- 
ference between the pair is very small), r(i) is not much 
different from ro and their difference essentially repre- 
sents random measurement errors. But because of the 
square, these errors do not cancel but add up after aver- 
aging over different pairs. 

So far we have shown that our measured pair disper- 
sions exhibit exponential growth in the dissipative range 
and power-law growth in the inertial range. However, 
these are manifested in different quantities, i.e. in (r 2 (t)) 
and ([r(i) - r ] 2 ) respectively. But in fact the pictures 
are consistent and there exists a crossover between the 
two regimes in both spatial and temporal scales. We 
note that Batchelor first discussed mean square separa- 
tion by using (r 2 (t)) -r§ rather than (|r(f) -ro| 2 ) [4|. In 
Fig. E£b) we plot several data sets with ro = 0.9, 1.3, 1.7, 
2.2 and 4. 377 using the original definition for pair dis- 
persion (r 2 (t)) -Tq [again normalized by /(r )t 2 ]. It is 
seen that the height of the these curves shifted down- 
ward systematically with increasing ro. As (r 2 (t)) = 
r 2 exp[0A2(t/r v )] for small values of t and ro, we plot 
(r 2 exp[0A2(t/T n )] -r 2 )/f(r )t 2 as the dashed red line 
in the same figure. It is seen that even in the high- 
resolution compensated plot the symbols agree excel- 
lently with the line. Note that the Taylor expansion 
of r 2 exp[0A2(t/T ri )] - r 2 with respect to time is domi- 
nated by 0A2r 2 t/T v for t/T v < 1. This can explain why 
the curves for small ro tilted up in the dissipative time 
range for mean-square separations compensated by t 2 . 
Figure [4jb) thus demonstrates the crossover from the ex- 
ponential to the Batchelor regimes both spatially (when 
ro varies from the dissipative to the inertial range of 
scales for fixed t <t T) ) and temporally (when t varies from 
smaller than r n to greater than r n for a fixed ro < 77). 

To summarize, we have made the first experimen- 
tal study of particle pair dispersions in buoyancy-driven 



thermal turbulence. In the dissipative subrange of scales, 
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FIG. 4: (Color online) (a) The mean square separation 
((r-ro) 2 } compensated by f{ro)t 2 as a function of normal- 
ized time t/T v for different initial separations measured at 
Ra = 1.3 x 10 10 {R\ = 84). From top to bottom, r =0.9r7 
(A), 1.3 77 (-4), 1.7 77 (►), 2.2rj (■), 4.377 (♦). There are 
18 data sets with ro = 6.5r7 ~ 65rj that collapse onto each 
other, which are represented by the solid circles . Inset: The 
black circles represent the value of ((r- ro) 2 )/i 2 at t = 2t v 
for different ro and the blue solid line represents /(ro) de- 
termined from Eulerian structure functions. Both quantities 
are normalized by er v . The red and green solid lines are 
/(r ) = fr 2 ,/^ 2 and/fro) = ±±Cr 2 ' 3 /rj 2/3 (C = 1.56 is the Kol- 
mogorov constant [lQ].), which are the dissipative and inertial 
range scaling predictions for /(ro), respectively. The purple 
dashed line shows the large ro limit of /(ro) = 6R\/\/T5. (b) 
[(r 2 (t)} - r§]//(ro)f 2 vs. t/T v . The open symbols here corre- 
spond to the closed ones in (a). The red dashed line shows 
(r 2 ,exp[0A2(t/T v )]-r 2 )/f(r )t 2 with r = O.977. 



our results show for the first time experimentally the 
existence of an exponential growth regime for the pair 
separation (r 2 (/;)}, which also yield the Batchelor con- 
stant B = 0.23 ± 0.07. For time t smaller and larger 
than to[= (ro/e) 1 / 3 ], respectively, the Batchelor and the 
Richardson-Obukhov scalings are observed in the mea- 
sured (|r(t)-ro| 2 ). The measured value of the Richardson 
constant is g = 0.10 ± 0.07. 
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Xu for helpful discussions. And gratefully acknowledge 
support of this work by the Research Grants Council of 
Hong Kong under grant CUHK404409. 
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